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I. INTRODUCTION 


Coherent-state path integrals for spin and for linear position and momentum degrees 
of freedom (and related phase-space path integrals) have been the subject of much study 
for over three decades now [l-9|, both for intrinsic reasons, and for their semiclassical limit, 
where they hnd application to many practical problems. Their mathematical subtleties have, 
however, prevented their widespread use, in contrast to the Feynman-position space path 
integral. For example, where in the Feynman integral the paths must be continuous but 
need not be differentiable, in the coherent-state case the paths need not even be continuous. 

In more recent years, steady pro gress has been made in understanding the semiclassical 


limit of such path integrals 


lOl, 


ll| . and in Ref. Il2| it was shown that when the so-called 


Solari-Kochetov (SK) correction is included, the resulting propagator in the spin case is 
free of j versus j + ^ arbitrariness, has the correct short-time behavior to O(T^), and is 
consistent under composition of successive propagators. In Refs. [iSl, [ij, this work was 
extended to coherent-state propagators for many particles and many spins. 


Insights gained from the above wor 


Sommerfeld quantization rule for spin 


lave led to the successful solution for the Bohr- 


15 


16 [|. an extension to the instanton calculus jlTj . 


and to a quantitatively correct explanation 


magnetic molecule Fe8(tacn)6 


18l. Il9| of the spin tunneling spectrum of the 


20| . Still, the answers for the propagators are obtained only 


by a careful examination of the discrete path integral, and casual application of meth¬ 
ods developed for the continuous-time Feynman path integral is fraught with errors. A 
continuous-time approach was adopted in Ref. 3, where it was found the path integral for 
the fluctuation determinant suffers from a global anomaly. The resolution of this problem 
again requires a careful examination of the discrete path integral, and it is shown that the 
anomaly is absent in a special gauge, whereby the Solari-Kochetov correction is automati¬ 
cally included. 

While these successes mean that the coherent-state path integral is no longer the hef- 
falump it once was, there is still some ambiguity in its conception. In particular, while it has 
long been known that the symbol (or c-number function) that plays the role of the Hamilto¬ 
nian in the classical action is not unique j^, how this nonuniqueness plays out in the hnal 
answer for the semiclassical propagator has not been properly explored. It is the purpose 
of this paper to do so, and in the process elucidate the nature of the SK correction further. 
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We will show that the SK correction is different, depending on the particular Hamiltonian 
symbol employed, but that the final answer is independent of this choice. Further, the final 
answer is best written using the Weyl symbol. As will become clear, this means that the 
formal, continuous-time coherent-state path integral is not only formal, it is also ambiguous. 
To give it meaning, one must return to the discrete path integral every time. 

Hints that the difficulties of coherent-state path integrals could be related to symbol- 


choice ambigui 


Refs. 


10 


11 


15 


ies (or what is the same thing, operator ordering ambiguities) may be seen in 


16l . Further, Kochetov 22| and Pletyukhov 23| noted that the SK correction 


could be written as the difference between the Q symbol and the Weyl symbol for the 
Hamiltonians, so that if, in contrast to previous papers which had employed the Q symbol, 
one employed the Weyl symbol in constructing the classical action, there would be no SK 
correction. Pletyukhov showed this result for a system with position and momentum degrees 
of freedom in generality, and for spin degrees of freedom within the Holstein-Primakoff 
approximation. The absence of the SK correction when one employs the Weyl symbol for 
the Hamiltonian suggests at first that it is simply clumsy to have worked with the Q symbol, 
and that if one uses the Weyl symbol from the outset, the correction will simply not arise 
in the first place. If true, this would be a nontrivial result since, as shown in Ref. Il2|, the 
correction arises from a global anomaly in the fluctuation determinant, and it is not clear 
how a change in the way the extremal action is expressed affects the fluctuations. Indeed, it 
is not clear how one would do the calculation with a general symbol in the first place. With 
this in mind, we calculate the propagator for particles in the P representation, following 


closely the derivation based on the Q representation in Sec. H of Ref. 


13 


Although the 


resulting change in the discrete path integral is seemingly minor, it leads to a nontrivial 
change in the final answer, and the SK correction now appears with the opposite sign. We 
then show that both this answer and the one from the Q representation are equivalent to 
each other, and to that for Weyl representation. We also show the analogous result for 
the spin case. In this we corroborate Pletyukhov, but we do not limit ourselves to the 
Holstein-Primakoff approximation, so our proof is completely general. The Weyl symbol 


for operators based on position and momentum degrees of freedom is of course classic 
but an analogous one exists for spin degrees of freedom too 
known. 


25 




29l |. although it is less well 


The plan of the paper is as follows. We present the results for the Weyl-representation 
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propagators for both particles and spin in the next section. This section also serves to 
introduce our notation, and to dehne principal terms. The P-representation calculation 
is done in Sec. HU and the equivalence of the Q-, P-, and Weyl-symbol-based answers is 
shown in Sec. IIVI The propagator for spin in the Weyl representation is derived in Sec. El 
In Secs. ED and IVIIl we turn to an examination of the continuum and discrete fluctuation 
determinants and their connection with the corresponding Jacobi equations with the goal of 
shedding more light on the global anomaly. Finally, in Sec. IVIIII we consider what happens 
when we try and evaluate the propagator for particles by working directly with the discrete 
action using the Weyl representation. Some essential facts about the Weyl representation 
(for both particles and spin) are collected in Appendix El 


II. PRINCIPAL RESULTS 


A. Propagator for particles 


For a particle with linear momentum p, coordinate q, and arbitrary Hamiltonian T-L, the 
propagator is dehned as 

K{zf,z,;T) = {zf\e-^^^\zi). ( 2 . 1 ) 

We have introduced here (unnormalized) harmonic-oscillator-based coherent states, 

|^) = e^“>), (^| = (0|e^ (2.2) 

with |0) and (0| being the normalized ket and bra for the ground state, and a and being the 
annihilation and creation operators. In Eq. (12.11) . Zi and Zf are arbitrary complex numbers. 
The Weyl form of the semiclassical approximation to the propagator K is 


/ 02 \ 1/2 

K^{zf, Zi] T) = h p- a; T) I exp [iS^{zf, zf, T)] 


dzfdzi 

Here, the classical action is given by 

iS^{zf, Zi] ^) = / + ^(0)^;*] -F 


-T r 


zz — zz 


[z{t),z{t)] 


dt, 


(2,3) 


(2.4) 


with H^{z, z) being the Weyl symbol for the Hamiltonian, and z{t), z{t) being the solution 
to the classical equations of motion 

dH^ 


dz 

dt 

dz 


dt * dz 


dH^ 


dt ^ dz ' 


(2.5) 

( 2 . 6 ) 
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with the boundary conditions z(0) = Zj, z{T) = Zf. 

We discuss the Weyl symbol at greater length in Appendix For now it suffices 
to recall the common textbook dehnition; is the c-number function obtained by sym¬ 
metrizing H in a and and then replacing these operators by the c-numbers ^ and 0 
respectively. 

The central point of the result fl2.3p is that it has no Solari-Kochetov correction. For 
comparison, when we use the Q symbol, the semiclassical approximation to K takes the 
form 12, 131 


K^{zf,z,-T)= z 


d‘^S^ 

dzfdzi 


1/2 


exp 




(2.7) 


The action is given by Eqs. (I2.4p - (l2.6p with the superscript W replaced by Q everywhere, 
with, additionally, the Q representation of the Hamiltonian 3^, being dehned by 

{zj+i\n\zj) 




Lastly, 


A® = 


{Zj+l\Zj) 

O^hQ 


( 2 . 8 ) 


(2.9) 


dzdz ’ 

and it is the term containing which we call the SK correction in Eq. fl2.7p . 

For completeness, we also give the answer for K when we employ , the P symbol for 
the Hamiltonian. We show in Sec. m that 

„• rT 


K^{zf,z,-,T)=[t 


Q2gP \ 1/2 


dzfdzi 


exp 


iS^{zf,Zi]T) - 


A^{t) dt 


( 2 . 10 ) 


All quantites here are the same as in Eqs. fl2.7p and 02.91) with the superscript Q replaced 
by P, and defined via 


n = 


d^z 


TT 


e-~^^H^{z,z)\z){z\ 


( 2 . 11 ) 


The signihcant point is that the SK correction enters Eq. 02.101) with a sign opposite to that 
in Eq. 02.7p . 

We will show in Sec. IIVI that , and are all equal up to the leading two terms 

in an expansion in h. 
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B. Propagator for spin 


The coherent-state propagator for a spin of magnitude j is dehned in parallel with that 
for particles: 

K{zf,zf,T) = {zf\e-^^^\z,). ( 2 . 12 ) 

The Hamiltonian is an arbitrary polynomial in the usual spin operators Jx, Jy and Jz- It 
follows that + Jy + is a constant of motion, which equals j{j -|- 1) for spin j. 

Further, \zi) and {zf \ are spin coherent states, dehned by 

k) = e^^-|j,j), (f| = (j,j|e^-^+, (2.13) 


with \iii) being the eigenstate of Jz with eigenvalue j, and J± = ± iJy. Again, the 

quantities Zi and Zf are arbitrary complex numbers which give the stereographic coordinates 
of the maximal spin projection direction in space. 

The Weyl-symbol-based semiclassical approximation to K is 

(2,14) 

Here, 

J = J + \- (2.15) 


This is reminiscent of the oft-stated prescription for the “classical” magnitude of the spin. 
Second, 


iS^{zf,Zi-,T) = + Zfz{T)) [l + z{0)zi)] + 


dt 


zz — zz 


J- 


1 + zz 


iH 


w. 


z, z 


( 2 . 16 ) 


Third, the path (z(t), z(t)J that appears in the action is the solution to the classical equations 
of motion, 

n -U rr'l2 I'T ^ f/W 

(2.17) 


dt * 2j dz 


(h _ . {l + zzfdH^ 

dt ^ 2? dz 


subject to the boundary conditions 2:(0) = Zi, z{T) = Zf. 

As the notation suggests, is the Weyl symbol for the Hamiltonian in the above ex¬ 
pressions. In contrast to the particle case, there is no simple analogue of the symmetrization 
rule for obtaining . Rather, it is dehned by the demands that the map from a spin opera¬ 
tor F to its Weyl symbol {z, z) be linear, covariant under rotations, yield a real c-number 
function for Hermitian operators, and, most importantly, obey the traciality condition 

d^z 


+ 71 j (l + ^;^)2 




(2.18) 
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for any two spin operators F and G and their corresponding Weyl symbols ^ . See Ref. 


29 


and references therein for details. A brief catalog of the resnlts most relevant to this paper 
is given in Appendix IA21 

Again, the signihcant point is that the form fl2.14|) of the propagator needs no Solari- 
Kochetov correction. By contrast, the answer based on the Q representation is 

1/2 


K^{zf,z,-T) = 


,{l + Zfz{T)){l + z{Q)zi) 525 'Q 
2j dzfdzi 


exp 


iS^{zf,Zi-,T) + - / A^{t)dt 

^ ./n 


(2.19) 


with 


iS^{zf,Zi-,T) = jln[(l + z/z(T))(l + 2 ;( 0 )zi)] + 


_ 

'.zz-zz Q ' 

1 ^ , iH^iz.z) 

Jo 

Gt 

+ 

—1 


= - 


1 f d {1 + zz)"^ dH^ d {1 + zzY dH^\ 


dz 2j dz ^ dz 2j 


dz ) 


( 2 . 20 ) 

( 2 . 21 ) 


and 

( 2 . 22 ) 

The qnantity A is the integrand of the Solari-Kochetov term, and is the Q symbol for 
Ti. The classical path obeys Eq. (I2.17p with in lien of . 

We do not bother to write explicitly; it would be completely parallel to Eq. (I2.19p . 
with the sign of the SK term reversed. This follows from what we do in Sec. IIIII and its 
extension to spin as indicated in that section. 


III. COHERENT-STATE PROPAGATORS FOR PARTICLES 


In this section, we consider the coherent-state propagator for particles. We show in 
Sec. IIII Al that the path integral for the propagator is not unique, and illustrate this by 
giving three different expressions for it. The hrst two are based on the Q and P symbols 
for the Hamiltonian, while the third uses alternating Q and P symbols. The semiclassical 


propagator studied in Refs. 




12h14| is the one based on the Q-symbol expression, and 


it contains the original SK correction. We will calculate the P-symbol-based propagator in 
Sec. IIII Bl where it will be seen that the SK correction arises with the opposite sign from that 
when the Q symbol is employed. The calculation starting from the mixed P-Q expression 
will be given in Sec. IVIII Bl 
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A. Setting up the path integral 


We expect on general grounds that, in any semiclassical approximation, K ~ exp(iS'), 
where S is the action for the classical path running from the initial state to the hnal state. 
However, this degree of approximation (analogous to the eikonal approximation in the WKB 
method) is too crude, for it ignores the conservation of probability. For that, one must 
include the next term in an expansion in powers of h. This correction generally takes the 
form of a pre-exponential factor ~ (d’^S/dzfdziY^'^, but of course we must hnd it more 
precisely. It is clear, however, that one must also calculate the exponent or eikonal correct 
to the leading two orders in an expansion in h. All these points are well known, but 
we dwell on them because the next-to-leading-order term is the source of all the trouble 
in all coherent-state-based semiclassical propagators and of Solari-Kochetov corrections in 
particular. Finding this term correctly is important as it is the one that assures conservation 
of probability. 

To calculate K, we divide the interval T into M slices of width A each: 


A = T/M, 


(3.1) 


where M S> 1, so that A is inhnitesimal and an expansion in A is permissible. We then 
write 

^ ... ^-mA factors). (3.2) 

Next, we insert a resolution of unity between every pair of adjacent factors in Eq. fl3.2p . 
Correct to order A, the propagator for one time slice is now evaluated as 


{zj+i\e "'^^\zj) = {zj+i\zj) exp[-iAH^{zj+i, Zj) + 0{A‘^)), (3.3) 


where 


3Q| 


H^{zj+i,Zj) = 


{zj+i\n\zj) 


Inserting the explicit expressions for overlaps such as {zj+i\zj), we obtain 


(3,4) 


rM-1 


K{Zf,Zi]T) 


jj- I (Pzj 


i=i 


TT 


exp 


‘SD. 


(3.5) 
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where the discrete action, is 


~ {zmZM-1 — Zm-1ZM-i) + {zm-IZM-2 — ZM- 2 ZM- 2 ) + ' ' ' 

M-1 

+ {Z2Z1 - ZiZi) + ZiZq - iA E H^{zj+i,Zj). (3.6) 

j=o 

Here, zm = %, and zo = Zi^ and it shonld be observed that 5'®^^ does not depend on zm and 
Zq for the simple reason that no snch variables have been dehned in the hrst place. 

We obtain a different expression for K based on the P symbol, if, again correct to order 
A, we write the jth factor from the right in the string fl3.2p as 


-,-iHA 


(fzi 


TT 


'^^—ZjZj—iAH^(zj, Zj ) 


Carrying ont this snbstitntion, we obtain 


K{zf,Zi;T) 


r M 

n 

■i=i ■ 


(Pz-i 


71 


^>\Zj){Zj\. 


exp 


where is another discrete action, given by 


(3.7) 


(3.8) 


'^^disc — {zm+iZM — ZmZm) + {zmZm-1 — Zm-1ZM-i) + ' ' ' 

M 

+ {z2Zi — ziZi) + Azq ~ iA ^ ^ 

i=i 


[Zj,Zj). 


(3.9) 


Again, zm+i ^ Zf, zq = Zi, and variables zm+i and zq do not exist, never having been 


defined. 

Equation fl3.9l) differs from Eq. fl3.6l) in two ways. The hrst is that we now have M 
integrations instead of M — 1. Since we are eventually going to let M —)■ oo, this change 
is insignihcant. The second difference is that is evaluated at Zj and Zj in the jth slice, 
whereas is evaluated at Zj+i and Zj. When we evaluate the extremal value of the action, 
we do so on a path where Zj+i — Zj = 0(A), so the second change would also appear to be 
inconsequential. Yet it is on precisely this difference that everything will pend, for it affects 
the essential properties of the two huctuation operators vis-a-vis their self-adjointness, or 
lack thereof. 

We obtain yet another expression for K if, instead of using all P’s or all Q’s, we alternate 
between the two. Let us consider the hrst two time steps starting with the state l^o) (^^o = Zi). 
We approximate propagation in the hrst step via i.e., we write 

^ j (3.10) 
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We propagate across the next time step by evolving the integrated-over state |zi) which 
appears above via , i.e., we write 

'■ d?Z2 _5 




ki) = 


--2^2-22 


\z2){z2\e 


vr 

<PZ2 


-2^2-22 


TT 


k2) 


g222:i g-*Aiy‘3 (22 ,21) 


(3.11) 


These two steps generate the following part of iSdisc'- 

^22^1 - + ^12^0 - 'i^{H^{z2, Zi) + H^{zi, (3.12) 

We continue in this way, alternating and H^. The resulting discrete action is 

**S'disc = {zm+iZM — ZmZm) + (zmZM-I — Zm-IZM-i) + ■ ■ ■ 

-^(^ 22:1 - ^ 12 : 1 ) + ^ 12:0 - ^ + H^{zj,Zj)). (3.13) 

The superscript A stands for ‘alternating’. 

It is clear that we can use and in any order, and thus obtain inhnitely many 
discrete path-integral expressions for K. We could also try and write {z'\e~'^'^^\z) in terms 
of the Weyl symbol using Eq. flA12D . extending the set of expressions even more. This 
immediately raises the question of how these different expressions will lead to the same 
semiclassical answer for K. We will address this question for the P representation in the 
next subsection, for the mixed P-Q representation in Sec. lVIIf^ and for the direct replaceent 
via the Weyl symbol in Sec. IVIII A[ First, however, let us see what happens if we take the 
formal continuous-time limit (A —)■ 0), and write the propagator as the path integral 

(3.14) 


with 


*5'fci = \ {zfz{T) + z{0)zi) 


dt 


zz — zz 


iH{z, z) 


(3.16) 


Not only is this equation merely formal, it is also meaningless, because H could stand for 
HP, H^, or something else, depending on which discrete path integral one starts with. This 


le corresponding path integral. If 
we will have to hrst specify what 


lack of meaning explains why there is an anomaly in t 
we try and work with Eq. 03.151) as was done in Ref. [12 
H{z, z) means, and, depending on that, the prescription for regulating the global anomaly 
will be different. This prescription will have to be obtained by examining the discrete path 
integral once again, so it seems that one is best off by working with the discrete form all the 


way, and eschewing the formal continuous-time form altogether 


32|. 
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B. P-representation propagator by integration by successive time slices 


In this section, we find the particle-case propagator starting with Eq. fl3.9p . We will do 
this nsing the method of Ref. Il3| since this method can be generalized to arbitrarily many 
particles and to arbitrarily many spins [l^. Since these references show how the extension 
to more than one particle or one spin is performed, we will show the calculation for one 
particle only, and leave the obvious generalization to many particles and many spins to the 
reader. We use the same notation, and focus on the changes that arise, so readers may wish 


to have a copy of Ref. 3 handy as they read along. 

The first step is to find the “classical” or extremizing path. The equations for this are 
essentially the same, and formally identical when we pass to the A —)■ 0 limit. The next 
step is to expand the action to second order in fluctuations around the extremizing path. 
Denoting the deviations in Zj and Zj from the classical path by rjj and fjj, the second variation 
of the action is 


r M 


r 2 qP _ 

^ ^^disc — 


2 ! 


Lj = l 


d 


d 




qP 


disc* 


(3.16) 


In terms of this quadratic form, the reduced propagator (the quantity multiplying the ex¬ 
ponential of the classical action times i) is given by 


K^ed{zf,Zi;T) = 


M 

n 

^■=1' 


d\ 


TT 


exp 


(3.17) 


As in Ref. Il3|, most of the derivatives in are zero. The exact expressions for the 

nonzero coefficients are slightly different, and they are now given by 


^jj ~ 




disc 






dz] 


dz] 


dz] 




a2 nP 

Djj = -i = 1 + iA 




dzjdzj 


f)2 oP 

Djj = = 1 + zA 


dzjdzj 




11 — z 


dzjdzj 

dzjdzj 


H {zj,Zj), 




[Zj,Zj), 


dzj^idzj 


= -l, 


'^disc 


dzidz. 


= - 1 . 


(3.18) 

(3.19) 

(3.20) 

(3.21) 

(3.22) 

(3.23) 


■j+i 
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Of these the hrst two are essentially the same as before (i.e., for the Q representation), but 
the last four are different. Clearly, Djj = Djj and Dj^j = 

The procedure at this point is to carry out the integrals time slice by successive time 
slice, and step three is to isolate the quantities that appear in the integral at the jth slice. 
We write this integral as 


/ 


(Prij 


exp 




n) y, 


(3.24) 


just as Eq. (2.34) in Ref. |l3|. The quantities rjj and r]j are the deviations in Zj and Zj from 


the classical path, 14- and 14- are row and column vectors given by 


^i = ~2 


y^ = 



0 



= 2^3+1 i 1 0 


V3+1 


(3.25) 


(3.26) 


G,= 


(3.27) 


and Gj is a 2 X 2 matrix that will be determined recursively. To avoid confusion with the 
time-slice labels, we label its elements with the letters “u” and “d” (for “up” and “down”), 
thus: 

Gj^uu Gj u^ 

Gj^du Gjdd 

The fourth step is to shift rjj and fjj so as to complete the square, and perform the 
integration for the jth slice. The shifts are given by 


, (3.28) 

7 ' . (3.29) 



[ij 7j) = Vj+i (1 0) G 


This leads, as before, to the consistency condition 


C — C 

^j,uu ^'i 


j,dd- 


This condition holds for j = 1, since 


ID-n aA 

\Dii Dll j 


(3.30) 


(3.31) 
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and Dll = -Dii- We shall see from the recursion found below that it holds for all j. The 
integral gives an overall factor of (detGj)“^/^, and a residual term in the exponent from 
completing the square, 


1 

2 








)ud' 


(3.32) 


Step hve is to examine the recursion relation relation for Gj and det(Gj). Equation (I3.32p 
implies that 

-^7+17 + 1 -^7+1 7+l\ (^7 


Gj + l — 


T+ii+1 j+ij+i 


) - (: 


)ud 


(3.33) 


This shows, hrst, that the consistency condition fl3.30p holds for all j. Second, as in Ref. 3, 
there is no meaningful recursion relation for det Gj, but there is one for the ud element Gud- 


Since {G^ )ud = —Gj^ud! det(Gj), this recursion relation is 


Gj+i^ud — Djj^i^j^i + (det Gj) Gj^udi 


(3.34) 


which, a priori, looks different from that in Ref. [iSl To see its explicit form, we note that 

(3.35) 


det Gj = 1 + 2iA 

OZiOZ 






which along with the expression for leads to 

a2rjP a2rjP \ 

This explicit form is the same as when we use the Q representation except that 
H^{zj,Zj) we have H'^{zj^i, Zj). Writing jA = t, and taking the limit A —>■ 0, 
turns into the Riccati differential equation. 


G 


j-\-l,ud 


= iA 


dz^.i 


(3.36) 

instead of 
Eq. (EJEI) 


~ iGud — B — 2AGud + BG\^, 


(3.37) 


which must be solved with the initial condition G'^ld(0) = 0. Here, 


A = 




B = 




B = 




dzdz ’ dz"^ ’ dz"^ 

The solution to this differential equation is, from Ref. Il3|, 

1 




(3.38) 


(3.39) 
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Here, 


6z{t) 

“W = m’ 


(3.40) 


which is a Jacobi held that describes how the classical trajectory for z(t) changes upon a 
change in the initial value of ^(0) while holding ;2(0) hxed. In particular, 


v(T) = i 


a'- 


dzfdz, 




-1 


(3.41) 


The quantity of greater interest to us, however, is not Gudif) but det G(t) (or detGj), 
since it is this determinant that we pick up from the integration at each time slice. The 
reduced propagator is 


M 




1 - 1/2 


(3.42) 


i=i 


Taking logs converts the product into a sum, which turns into an integral in the limit A —>■ 0. 
We found det Gj in Eq. (I3.35p . Hence, 

„• rT 




dt [2AP{t) - B{t)GUt)], 


(3.43) 


which differs from Ref. Il3| in the extra hrst term, 2A{t). Feeding in the solution (I3.39p . we 
obtain 


= -2 


A^{t) dt — - Inn(T). 


(3.44) 


Hence, the hnal answer for the propagator in the semiclassical approximation is, as advertised 
before. 


K^{zf,z,-,T)= it 


dzfdzi 


1/2 


exp 


iS^{zf,Zi]T) - 


A^{t) dt 


(3.45) 


IV. EQUIVALENCE OF THE PARTICLE PROPAGATOR IN DIFFERENT REP¬ 
RESENTATIONS 


Our goal in this section is to show that Eqs. fl2.10p and fl2.7p are equivalent, and to write 
the propagator using the Weyl representation. 

If we look at Eqs. flAlbp and flA20p it seems that we can replace and by 
in and and delete the Solari-Kochetov correction. This will turn out to be correct, 
but there is one subtlety which we must hrst mind. The path (^z{t),z{t)) which appears 
in the action is obtained by solving the equations of motion fl2.5p and fl2.6p but with 
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instead of . Let us temporarily denote the path by z^{t)) and the action by 

S^[z^(t),z^(t)] to emphasize this fact. Let us likewise denote the classical path based on 


by [z^{t), z^{t)y Since 

H^{z,z) = H^{z,z)x {l + 0{h)), (4.1) 

it follows that 

Z^{t) = z^{t) X (1 + 0(h)), z^{t) = z^{t) X (1 + 0(h)). (4.2) 

The action is, however, an extremal value. A small change in the path therefore changes 
the action only in second order. That is 

S'’[5"'(«),z“'(i)] X (1 + 0(S^)). (4.3) 

By the same argument, 

S^[z^{t),z^{t)] = S^[z^{t),z'^{t)] X (l + 0(h2)), (4.4) 

= S^[z^{t), z^{t)] X (l + 0(h^)). (4.5) 


Since our goal in calculating the semiclassical propagator is to obtain it correctly up to the 
hrst term in relative order h, these changes are beyond the accuracy to which we are working, 
and may be neglected. They may be similarly neglected in the prefactor (d'^S^/dzfdziY^"^. 
To this order of accuracy, therefore, . By the same argument, . 


V. SPIN PROPAGATOR IN THE WEYL REPRESENTATION 


In this section we shall give the propagator for spin following the ideas developed in the 


previous sections for particles, and notation developed in Ref. 


29 


Our aim is to obtain 


Eq. (I2.14p . starting with the previously obtained result, Eq. (I2.19p . i.e., to rewrite in 
terms of the Weyl symbol H^{z, z) in parallel with Sec. IIVI As a preliminary step, we hrst 
discuss the stereographic variables 2 ; and z, which are often more convenient descriptors of 
the phase-space sphere than the orientation h. If the spherical polar coordinates of h are 
taken as then 

(5.1) 


2 ; = tan 


2 ; = tan |e 


15 







The spin coherent state |h), which is the state with maximum spin projection along n, i.e., 

J ■ n|n) = jjh), (5.2) 

can clearly be obtained from the state with maximum projection along z, i.e., by 

applying a rotation. When the requisite rotation operator is written in terms of 2 ; and z, its 
action on \j,j) can be cast in the form fl2.13p up to a multiplicative constant, i.e., 

|n) oc \z). (5.3) 

This result is a proportionality rather than an equality because, as dehned in Eq. fl2.13l) . the 
state \z) and its dual bra {z\ are not normalized; rather 

{z\z') = {I + zz'fK (5.4) 


The resolution of unity now takes the form 


2j + 1 r <fz 

71 J {1 + zzy^+'^ 


(5.5) 


The beneht of using unnormalized states and z, z variables is the same as for particle 
coherent states: Matrix elements are analytic in and 0 , and we can exploit analyticity to 
simplify many calculations. 

Next, we note that, by Eq. flA30p . 

H^{z,z) = + (5.6) 


where £ = —z(h x Vn) is the angular momentum operator on phase space (not the Hilbert 
space of the states \j,m)). In terms of z and z [33|, 


= -(1 + ^; 2 ) 


dzdz' 


(5.7) 


Hence, Eq. flA30p may be written as 

\Z,Z) ti yz,z) 

which is correct up to relative order 1/j. The next step is to write 


(5.8) 


A{t) — Aiit) + A2(t), 


(5.9) 
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where 


= 

A2 = 


(1 + zzf 

2j dzdz ’ 

(l + zz) f dH^ _dHQ\ 


2j 


dz 


+ z- 


(5.10) 

(5.11) 


dz J ' 

In Ai, we may replace 2j by 2J in the denominator since A is already of order 1/j relative 
to S, and we do not care abont errors of relative order 1/j^. Thus, 


Mt) = 

For A 2 , we recast it by using the equations of motion. Thus, 


(5.12) 


A2it) = -i 


zz — zz 
1 + zz ' 


(5.13) 


The terms Ai and A 2 can be combined, respectively, with the second and the hrst terms in 
the integral in Eq. fl2.20p to yield 


dt 


.zz — zz 
' 1 + zz 


— iH^{z, z) 


+ 


A{t) dt = 


dt 


dt 


zz — zz 




-nl + 7v)if^(^,^) 

1 + zz V 4j / 


.„Wf- X 

J—— - iWAz,z) 

1 + zz 


(6.14) 


Next, we observe that the coefficient of the explicit boundary term in the action fl2.20p can 
be changed from j to j by lifting the corresponding term in the prefactor into the exponent. 
In this way, we obtain 

• Q2gQ \ 1/2 


K^{zf,Zi-T) = exp{jln[{l + Zfz{T)){l + z{0)zi)]) 


X exp 


dt 


^zz — zz 


— -iH^{z,z 


1 + zz 


(5.15) 


We can now make two further changes which only affect our answer to relative order 1/j^. 
First, we can employ the same argument which led to Eq. (14. 3 p to replace the path used to 
calculate the action be the one based on instead of Second, we can replace the j 
in the prefactor by J. This gives us the Weyl-symbol-based propagator, Eq. fl2.14p . 
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VI. THE CONTINUUM FLUCTUATION OPERATOR AND THE CONTINUUM 


JACOBI EQUATION 


While the action in the formal continuum limit, S'fd, is ambiguous it is still useful to 
consider the reduced path integral, 




fcf = / [dz] [dz]exp{iS^S{4v,v]), 


where 


= “5 / (>)(*) vW) ®tei ^ I 


( 6 . 1 ) 


(6.2) 


d^icl — 


(6.3) 


'0 ^ ^ \r]{t) 

with Dfci being the fluctuation operator 

-idt + A(t) B(t) 

B(t) idt + A{t)^ 

acting on paths that obey the constraints t]{ 0) = fi{T) = 0. The ambiguity in the formal 
continuum limit shows up as follows. As found in Ref. 3, because of the global anomaly 
in the path integral the operator "Dfci has no eignefunctions, not even one. Thus, 

detPfci cannot be dehned as the product of the eigenvalues of Pfd. For the same reason, 
the standard method for hnding this determinant based on solving the associated Jacobi 
equation also fails. 

To explain the nature of this failure, we now describe the Jacobi-equation-based method. 
The classical equations based on the continuum action are 

(6.4) 

These equations have to be solved with the boundary conditions 2:(0) = Zi, z(T) = Zf. The 
other boundary values, z(0) and z(T) are not fixed, but emerge from the solution and may 
thus be regarded as functions of Zi, Zf, and T. If now we use the value of ^(0) so found and 
Zi to solve the classical equations of motion as an initial value problem, we will recover the 
classical solution for z{t) and z{t). If we change the initial values to Zi and ^(0) -1- e, where e 
is inhnitesimal, the solution to the initial value problem will deviate from the previous one 
by terms of order e in leading order. Denoting the deviations in z{t) and z{t) by eu{t) and 
ev{t) respectively, we hnd the Jacobi equations, 

-idt + A{t) B{t) \ lu 
Bit) idt + Ait) / \v 


dz .dH 

dt dz ’ 


dz .dH 
dt * dz 


= 0 , 


(6.5) 
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with initial conditions m( 0) = 0, n(0) = 1. (Since we will not need it, we do not bother 
writing the full Jacobi system allowing for variations in Zi also.) The Jacobi-based method 
says that 


f.Q, 


det'Dfci = v{T). 


( 6 . 6 ) 


As found in Ref. [12, Eq. (16.6p is incorrect and should be multiplied by an undetermined 
phase factor, . This phase factor is the SK correction, which we now know differs for the 
P and Q representations, while T>ic\ is superhcially the same in the two cases. 

For completeness, and to enable the reader to understand the answers fl2.3p . fl2.7p . and 
fl2.1Up for the propagator, we mention that 


v(T) = ( 




dzfdzi 


The proof is standard. See, e.g.. Sec. 4 of Ref. 


12 


or Ref. 


34 


(6.7) 


VII. THE DISCRETE FLUCTUATION OPERATOR AND THE DISCRETE JA¬ 
COBI EQUATION 

In this section, we return to the discrete path integral and, for both the P and Q rep¬ 
resentations, examine the fluctuation operator by writing it as a tridiagonal matrix. The 
same operator determines the Jacobi equation. We will show that unlike the continuum 
case, the determinant of the discrete operator is not simply equal to the solution to the 
discrete Jacobi equation but also contains an SK correction. We will see why the correction 
differs between the two cases (P and Q). We will further see that the solution to the discrete 
Jacobi equation tends to the continuum solution in the limit A —)■ 0. This shows the precise 
way in which the equality of the Jacobi held and the huctuation determinant breaks down 
in this limit. 

This method cannot be extended (at least we do not know how) to more than one particle 
or spin, but the insights it provides as described above still make it worth presenting. 


19 














A. The source of the SK correction 


For either Eq. fl3.9p or Eq. fl3.6p . we can write (taking the number of intermediate inte¬ 
grations as M in both cases) 

hi 
hi 
h2 
h2 


*<5^*S'disc[h>h] = - 


2 L 


hi hi h2 h2 


hM hM 




disc 


hM 

hM 


where "Pdisc is the discrete fluctuation operator (or matrix) 


7~) _ 

■^disc 


Ai 

-Dll Dll 

Di2 



D21 

D22 D22 

D22 D22 

b 

to 

COI 



D32 


Dm-im 



Dm M-1 

Dmm DjPh^ 

Dm M Dm m _ 


(7.1) 


(7.2) 


Note that we have reordered the r^’s and i^’s in the row vector, as this makes Pdisc a manifestly 
symmetric and tridiagonal matrix, albeit complex 35|. The reordering leads to M additional 
factors of —1 when the determinant is evaluated, so that 


/^red= [(-!)"" detPdisc] 


- 1/2 


(7.3) 


Next, let us examine the discrete Jacobi equation. The equations for the classical path 
that follow from the discrete action can be written as 

^(-*^disc) = 0, _(-*S'disc) = 0, (j = l,2,...,M). (7.4) 

To derive Jacobi equations from these, we would like to treat zq = Zi and zq as initial values. 
This, however, is meaningless as there is no such variable as Zq. Instead, we must take Zq 
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and zi as the initial valnes, the latter being regarded as determined by Zi and Zf = zm+i- 
We now keep Zi nnchanged, and let zi ^ zi + e. Let ns denote the changes indnced in Zj 
and Zj by 

6zj = euj + O(e^), 6zj = evj + 0{e^). (7.5) 

Performing the necessary variations, we obtain, 


E( 


d d 

Uk^ -h Vk^ 

OZk OZk 


dZjj 

-iOSf^iscj Ozj 


)-l 





(7.6) 


Most of the derivatives on the left vanish. When only the nonzero ones are kept, we obtain 


DjjUj + DjjVj + Dj^jVj+i = 0, (7.8) 


nsing the dehnitions of the varions D^s. [These dehnitions can be read off by considering 
only the hrst member of each of the eqnations fl3.18p - fl3.23p and deleting the ‘P’ snperscript.] 
These eqnations hold for j = 1,2,..., M, and we mnst take uq = 0, vi = 1. They then 
determine Uj for 1 < j < M, and Vj for 2 < j < M + 1, all of which are meaningfnl 
qnantities. We now observe that we can rewrite them in the form 


Dll Dll 





Vl 


0 

Dll Dll 

Di2 




Ui 


0 

D 21 

D 2 2 D 2 2 




V2 


0 


D 2 2 D 2 2 

D 23 



U 2 


0 


D 32 










Dm-im 







Dm M-i 

Dmm 


vm 


0 




Dm M Dm m _ 


. “at . 


_ Vm+1 _ 


We have nsed the initial condition on uq {uq = 0), bnt not on ni, leaving it as an arbitrary 
qnantity instead. In fact, the way this eqnation is written snggests that all Vj, Uj (j = 
1,..., M) are determined in terms of vm+i- By demanding that vm+i mnst be chosen snch 
that fi = 1, however, we once again obtain vm+i explicitly. 

The matrix that appears in Eq. fl7.9p is of conrse none other than T’disc- By Cramer’s 
rule, therefore. 


detC 
det fPdisc 


(7.10) 
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where, 


0 Du 

0 Hii 

Di2 



0 D21 

D22 D22 



0 

D22 D22 

7^2 3 



D32 






Dm-im 

0 


Dm M-i 

Dmm Dj^j^ 

Vm+1 



Dm M Dm m 


To evaluate detC, we expand it by the hrst column, obtaining. 


(7.11) 



Du 






Dll 

Di2 





D21 

D22 

D22 



det C = —vm+i det 


D22 

D22 

D23 




Ds2 


Dm-im 




Dm M -1 

Dmm Dmm 


(7.12) 


The matrix that remains is lower triangular, so its determinant is just the product of the 
diagonal entries. Anticipating future minus signs, we dehne the quantity 

M M 

Tsk = (-l)"'"TI^sn^J-iJ- 

i=i i=2 


in terms of which 


detC — (—Ij^UAT+ihsK 


(7.14) 


and 

Vi = (- 1 )^ , vm+i- (7.15) 

det Pdisc 

Setting ui = 1, and rearranging, we get 


(—l)^det T>disc — hsK'^M+i- 


(7.16) 


This is the correct discrete replacement of Eq. fl6.6p . The factor Tgx need not be unity, in 
which case we have a nonzero SK correction. 
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The final step is to take the A —)• 0 limit of Eq. fl7.16p . By definition, the left hand side 
turns into the continuum fluctuation determinant, and its inverse square root will give us 
Kred- It remains to see what happens to vm+i and Tsk on the right hand side. It is simplest 
to do this separately for the P and Q representations. Before turning to this, however, it 
pays to rewrite the general Jacobi equations, fl7.7p and fl7.8p . as the 2x2 matrix recursion 
relation. 




Ui 




^3J 


Uj-i 


(j = l,2,...,M), 


(7.17) 


y ^jj ^3+ij J J \ ^ / \ ^^ / 

with the initial conditions Mq = 0, Ui = 1. 

It is of course also possible to evaluate det Udisc by writing a recursion relation for suc¬ 
cessive diagonal sub determinants of Udisc (as may be done for any tridiagonal matrix). The 
SK corrections can then be obtained by examining the relationship of this recursion relation 
to Eq. fl7.17l) . We shall not follow this route. 


B. Application to P representation 

Let us consider the discrete Jacobi equation first. We feed the explicit values of the ZJ’s 
from Eqs. fl3.18p - fl3.23p into Eq. fl7.17p . and abbreviate 


A(zj,Zj) — Aj, B(zj,Zj) — Bj, B(zj,Zj) — Bj. 


(7.18) 


We find that 


-{l + iAAj) 0 


Un 


-1 i/S.Bj \ / Uj-i 


iABj —1 / \vj+i } \ 0 1 J- iAAj / \ 

Solving for the column vector on the left, and dropping terms of O(A^), we obtain 


1 — iAAj —iABj \ I Uj-i 


(7.19) 


/ Uj \ 



1 1 


iABj 1 J- iAA 


(7.20) 


3/ \ ^3 


It is immediately apparent that in the continuum limit (A —)■ 0, M —)■ cx), with AM = T 
fixed), this recursion will turn into the continuous-time Jacobi equation flO.Sp . Since the 
initial conditions are also identical, it follows that 


lim vm+i =v{T). 

A->0 


(7.21) 
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The second step is to evaluate Tsk- We have, 


~ 1 + iA.Aj, Dj-ij — 1 . 


(7.22) 


M 


Hence, 

Tsk = (-1)^"^ X JJ(1 + iAHj) x (-1)^"^ cs exp(iA^Hj). 

i=i j 

Collecting together the above results, we hnd that 

lim (—1)^ detPdisc = expfi [ A(t)dt)v(T), 

A-S-O V /n / 


so that 




- 1/2 


(7.23) 


(7.24) 


(7.25) 


The extra exponential factor is the SK correction, and we see that we have the correct sign 
for it. 

C. Application to Q representation 


We now repeat the previous subsection’s arguments for the Q representation. The relevant 
D coefficients are given in Eqs.(2.16)-(2.21) of Ref. Il3|, and we redisplay them here for ready 
reference: 


^jj ~ Djj — iABj, 

where now, 

Aj = A(zj, Bj = B{zj, Zj_i), Bj = R(%i, Zj). 

The discrete Jacobi equation now reads 

T 0 \ I Ui \ I—1 + iAAj iABj\ / Wj-i 


(7.26) 


0 


-1 


(7.27) 


(7.28) 


iABj —l + iAAj_^_il yuj+iy . w ^ ^ 

Again we solve for Uj and Vj+i to 0(A). In the process, we also replace the AA^+i term by 
AAj, since the difference is O(A^). In this way, we get 


Un 




1 — iAAj —iABj 
iABj 1 + iAAj 


Uj_i 


(7.29) 
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which is formally the same as in the P case. For the same reasons as given there, we again 
get 

\im vm +1 = v{T). (7.30) 

A^O 

Next, for Fsk, we have 


M-l 


Fsk = 

(-1)" 1 X (1)^ X JJ(-l + iAAj) ^exp(-iA^A^). 

i=i j 

(7.31) 

It follows that. 




lim (— l)^detDdisc = expf— z / A{t)dt)v{T), 

A^O ^ Jo ' 

(7.32) 

and 

. n'j' 



K?ed = exp A{t) dt^ [v{T)] 

(7.33) 


The extra exponential factor is the SK correction for K^. We draw the reader’s attention 
to the sign. 


VIII. DIRECT EVALUATION OF THE PARTICLE PROPAGATOR IN THE 
WEYL REPRESENTATION 

Our goal in this section is to try and evaluate the propagator using the Weyl represen¬ 
tation for "H from the very start. One way to try and do this is to write the infinitesimal 
time-evolution operator in terms of using the Weyl kernel W(f, z) as in Eq. flA5D . 

We shall see that this way does not work. The other way, which does work, is to alternate 
P and Q representations, building on Eq. (13.131) . 

A. Mapping via Weyl kernel 

First, let us use Eq. flA5l) to write the infinitesimal time-evolution operator as 



where by [X]-ws we mean the Weyl symbol of the operator X. Using Eq. flA12p . we may 
write the propagator for one time slice as 
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To find \e we expand e 

H^{z, z), and that for is 


-mA 


in powers of A. The symbols for 1 and 'H are 1 and 







4 

dz‘^ dz‘^ 

V dzdz / 


+ ■■■ . 


( 8 . 3 ) 


The important point here is that the correction is of relative order and not h. (We will 
show the powers of h relative to the leading term explicitly in this section.) Hence, when 
we reexponentiate the series, we hnd 


[e-^^^]^^=exp{-iAH^{z,z) + 0{h‘^A‘^)). ( 8 . 4 ) 

For the time-slice propagator, we get 

{z2\e-^^^\zi) = J ( 8 . 5 ) 

with 

<F(f, Z 2 , zi) = -iAH^(z, z) - 2{z - Z 2 ){z - zi) + O(h^A^), (8.6) 

The natnral procedure at this point is to evaluate the integral over z and z semiclassically, 
i.e., by steepest descents. Let us denote the critical (saddle) point by Zc and Zc, and partial 
derivatives by subscripts. Setting $2 = = 0, we hnd 

- -^AHY{zc,z,) + 0{n^/‘^A^), ( 8 . 7 ) 

Zc = Z2 -^AiLf {zc, Zc) + 0(h®/^A^), (8.8) 

where we continue to show powers of h explicitly. Hence, denoting the critical value of <F by 
<hc, we have 


= -tAH^{zc,Zc) + lA^hHY{zc,Zc)HY{zc,Zc) + 0{h^A^), ( 8 . 9 ) 

the additional error terms introduced at this step being of order h^A^. Since and Zi differ 
from Z 2 and Zi by terms of order Ah^^^, it is reasonable to perform a second expansion in 
A. When this is done, we hnd (the sign of the A^ term should be noted) 

= -tAH^{z2,zi) - \A^nHY{ z2,Zi)HY{Z 2,zr) + 0{A^h^). (8.10) 


The next step is to perform the Gaussian integral over the small deviations from the 
critical point. Dehning rj = z — Zc, f] = z — Zc, 'we have 


$ 


4 >, 


+ ^ + ‘^^zzm + ^zzt/^ h — , 


( 8 . 11 ) 
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with 


= -2-i^hHY,, ( 8 . 12 ) 

. 

At this point, it is better to leave the derivatives of evaluated at Zc, Zc- The integral 
gives us the inverse square root of the determinant of this quadratic form, which equals 

1 - + 0(A"A^). (8.13) 

Evaluating ifT at ^ 2 , zi incurs a further error of the same order, i.e., 

Putting all these pieces together, we hud, eventually, 

= exp (2221 - tA[H'^ + + 0(lW)) . (8.14) 

The order A terms combine to form H^{z 2 , Zi), so that we get, as before, 

(Z2|e-*^^|zi) ^ (8.15) 

but now we know the order of the terms omitted. 

We thus see that it does no good to start with , since that entails auxiliary integrations 
for each time slice, which when performed will lead to the same fluctuation determinant 
over the non-auxilliary variables as before. Only the Q and P representations allow us to 
dispense with auxilliary integration variables, but then the fluctuation determinant leads to 
SK corrections. 


B. Alternating P and Q representations 


We now start with the discrete action fl3.13p obtained by alternating P and Q representa¬ 
tions. This is an obvious thing to try, since when we determine the extreme or the classical 
value of the action, the alternating if^’s and iP'^’s will combine to produce as A —)■ 0. 
The expectation is that the reduced propagator will then be free of any SK correction. 


The only nontrivial part of the calculation is the in 
we do by successive time slices as in Sec. [TTIl and Ref. 


egration over the fluctuations, which 


13l. picking up factors of (detGj) 


at each step. Suppose the step from j to j -|- 1 is of type P, and the next step is of type Q. 
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Then, adopting a notation for the derivatives of and analogous to that in Djj, etc., 
the part of 6‘^{iS) which involves the fluctuations at these steps is 


= ■■■- - (1 + 

+ HfH + (1 - 

- ’fe+i5,+i + %+i’7j+2 + ■ ■ ■ ■ (8.16) 

We do not show the values of ^ and at which the derivatives of and are evaluated 
explicitly; they are {zj,Zj) and {zj+i,Zj), respectively. Since we will eventually let A —)■ 0, 
these arguments will take values on the classical path z(t), z(t) with t = jA. 

For the integration at step j (over rjj and fjj), we need the matrix Gj as well as the vectors 
Vj and Vj. We have 

G, = f ^ \ , ( 8 , 17 ) 

V, = 1(1 - , Vj = 1(1 - iAHf-^) (l o) . (8.18) 

The matrix element Gj^ud is of course unknown, having been modihed as a result of the 
previous integration steps. The integration at this step produces a residual term from 
completing the square equal to 

5(1 - W.(G-‘)„, = -5j^(l - + 0(A=), ( 8 , 19 ) 

and a determinantal factor (detGj)“^/^. It is apparent that 

det Gj = l + 2iAHl. - iA {Hf- + H^) Gj^^d- (8.20) 


The result of step j is that the matrix Gj+i equals 

^ , (8.21) 

0 1 J 

where 

= iAH^— + -^(1 - 2iAH^—)GTud- ( 8 . 22 ) 

J + i.na j+lj+1 Q. \ ],J+IZ ^ > 




Further, 


= ^Vj +2 , Vj = ^r]j +2 (l o) . (8.23) 

The integrations at step j + 1 thus produce a determinantal factor (det = 1, and 

a residual term from completing the square equal to 

2%+2(^j+i)ud ^ ~2^j+2^j+i,ud- (8.24) 

Hence, 


G, 


j-\-2,ud 


*^-^j+2j+2 


G 


j-\-l,ud 


3+2,j+2 j+lj+l 


+ 1 


2tAHf— - 2tAHP^ + tA{Hf^ + 


(8.25) 

(8.26) 


This is the recursion relation desired, since we have now integrated over the complete repeat 
pattern. If we let A —)■ 0, it will be seen that the symbols and always appear in the 
combination 

H^{z, z) + H^{z, z) = 2H^{z, z){l + (8.27) 

The differential equation for Gud is 


-iGud = B^ -2A'^Gud + B^G: 


w, 


^udi 


(8.28) 


where we have added a superscript W to show that the symbol for the Hamiltonian that 
is involved is . This differential equation is the same as before with the same initial 
conditions. Hence, 




(8.29) 


B^{t) V ' ' V. 

Finally, we need the product of all the determinants detGfc. Keeping in mind that the 
determinant from every other step is unity, we have 

. M 


In Kj-ed 


- - ^ln(det Gk 


1 

'2 

1 

2 


k=l 

•T. g2ffP 


'0 


L dzdz 


Gud{t) 


rT ^ 


A^{t) 


d2HW 

A^{t)+i^ 

V 


dt 

dt. 


(8.30) 
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Now, = 0{h), which may be neglected since the term we are discussing is already 

the hrst correction in powers of h. Hence Ini^^red = i.e., 

q2 qW \ 1/2 


— I ^ 


\dzfdzi) 


( 8 . 31 ) 


which has no SK correction. 
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Appendix A: Review of P, Q, and Weyl symbols for particles and spins 


1. Mapping for particles 


For a massive particle in one spatial dimension, the Q and P symbols of the Hamiltonian 
l-L are dehned by 

{z\n\z) 


H^{z,z) = 

H = 


{z\z) 

f —e-~^^H^{z,z)\z){z\. 
rc 


(Al) 

(A2) 


(Analogous dehnitions apply to other operators.) In Eq. flA2l) (fz is shorthand for dxdy, 
with X and y being the real and imaginary parts of z. The reason for the factor inside 
the integral for this equation and the {z\z) denominator of the previous one is that as dehned 
in Eq. fl2.2l) . the states \z) and {z\ are not normalized; instead 


{z\z') = 


(A3) 


The resolution of unity therefore takes the form 

1 = J ^e-~^^\z){z\, (A4) 

with the same extra factor. The advantage of using unnormalized states is that off- 
diagonal matrix elements such as (z|'H|z') can be obtained from the diagonal one, (zlHlz), 
by appealing to analyticity. In practical terms this means that we merely replace 2 ; with z'. 
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For the Weyl symbol, we follow Weyl himself 2^, and dehne 

H^{z,z) =Tt{HW{z,z)), 

where W(^, z) is an operator-valned kernel given by 

d?w 


}V(z,z) = 


(wa'^ —wa) —(wz—wz) 


7T 


(A6) 


(A6) 


This dehnition pnts the familiar symmetrization rnle for Weyl ordering of operators on a 
broader footing and can be shown to rednce to that for simple examples snch as a^(a^)^. The 
Weyl symbol for any other operator is dehned analogonsly. By letting w —)■ —w, w —)■ —w 
in Eq. (1A6D . we hnd that [W(f, z)]^ = W(f, z). Hermiticity of "H then implies that H^{z, z) 
is real. 

One possible inverse of the transform flA5l) is given by 

d?z 


n = 


TT 


H^(z,z)W{z,z). 


(A7) 


It is straightforward to show this resnlt by using the identities (themselves easily shown) 


TiW{z,z) = 1 , 

Tt (yV(z, z)yV(z', z')) = 7r6^'^\z — z'), 


(A8) 

(A9) 


where by 5^‘^\z) we mean (5(Re 2;)(5(Im ^). In fact, the inverse flA7p is unique, as may be 
shown by taking the trace in Eq. flASp in the complete set of position states. The matrix 


elements of the kernel W are not difficult to find, and the trace takes on the form of an 
ordinary Fourier integral, which may be inverted to obtain an expression for the position- 
space matrix elements of 'H in terms of the function . This expression is easily seen to 
be identical to the one implied by Eq. flA7p . Thus Eq. flA7p not only implies but is also 
implied by Eq. 
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Next, we recapitulate the relationship between the Q, P, and Weyl symbols 
the matrix element 

(^21771^1), (AlO) 

which is nothing but {z 2 \zi)H^{z 2 , Zi). (See Eq. fl3.4p .) Writing 77 in terms of H^{z,z), 
and {z 2 \zi) = we obtain 

<Pz 


HQ{z2,zi)=e-'^^^^ I —H^{z,z){z2\W{z,z)\zi). 
. vr 


(All) 
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Using Eq. 


, we get 


cPw 


(f 2 |>V(^, z)\zi) = / - exp{z 2 W - wzi - ^ww + Z 2 Z 1 ) 


TT 

cPw 

71 


exp(-|ww + w{z 2 - z) - w{zi - z) + Z 2 Z 1 ) 


= 2 exp(^- 2 {z 2 - z){zi - z) + Z 2 Z 1 ). 

We now feed this result into Eq. flAllD while at the same time dehning 


(A12) 


T] = Z - Zi, T] = Z - Z2. 


(A13) 


We thus get 

H^{z2, ^ 1 ) = 2 y ^e-^^m^{z2 + r],z, + 7]). (A14) 

If we now Taylor expand in powers of r] and fj, it is easy to perform the resulting 
Gaussian integrals. Retaining the hrst nonzero correction, we get 

H^{z2, z,) = H^{z2, Zi) + --^^H^{z2, Zi) + • • • . (A15) 

The second term in this expansion is in fact of order h relative to the hrst. One can see this 
point by writing the quantities z and z in terms of dimensionful position and momentum 
variables, and noting that ^ and ^ both contain a factor of gy transposing this term 

to the left hand side, and using the same equation recursively, we hnd that 

H^{z2, z^) = HQ{z2, z,) - 1-^H<^{z2, Z,) + ---. (A16) 

2 OZ 2 OZ 1 

To relate and , we substitute Eq. flA2p in Eq. flASp . and obtain 

H^{z2,zi) = T:t{HW{z2,zi)) 

= j —e-'^^H^{z,z)T:i{\z){z\W{z2,zi)) (A17) 

Now, 


i:i{\z){z\yV{z 2 ,Zi)) = {z\W{z 2 ,Zi)\z) 

= 2exp(^—2{z — Z 2 )(z — Zi) + zz^, (A18) 

where the last result is obtained from Eq. flA12p with the exchange {z 2 , zi) y-)- (7, z). Feeding 
it into Eq. flA17p along with the dehnitions flAlSp . we get 

—+ ■ 7 , 2^1 + ??)• (A19) 


H^{z2,z,) = 2 


71 
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We now Taylor expand in powers of rj and fj just as done above, and integrate over rj 
and fj. Again retaining only the first nonzero correction, we get 

1 (9^ 

{z2, Zi) = H^{z2, ^i) + - „— H^{z 2, Zi) + ■ ■ ■ . (A20) 

2 OZ2OZ1 

For the particle case, we can achieve a more general correspondence between operators 
and phase-space functions by extending the definition of the Weyl kernel to 

W^'^\z,z) 

where — 1 < a < 1. The mapping is then given by 

H^‘^\z,z) = Tr (A22) 

n = J ^H^^\z,z)W^-^\z,z). (A23) 

The cases of Q, P, and Weyl mappings correspond to a = 1, —1, and 0 respectively. We 
shall not employ this general definition, but shall work with the P and Q mappings in the 
form given earlier. 


d^w 


^{wa^ —wa) ^—(wz—wz) ^azz/2 


TT 


(A21) 


2. Mapping for spin 


Let us hrst discuss what we mean by the Hamiltonian of a spin system. For a particle of 
spin j, the most general Hamiltonian (or any other operator) can be written as 381 


2j ^ 

C£m3^rm(J); (A24) 

£=0 m =-£ 

where Qm are arbitrary c-number coefficients, and Wm are spherical harmonic tensor op¬ 
erators defined via the operator analogue of the Herglotz generating function for spherical 
harmonics. The classical phase space can be taken as a sphere of fixed radius (which may 
be taken as 1, or j, or j + ^, whichever is most convenient). As discussed in Sec. El the 
variables z and z that we employed in Sec. Ill Bl are stereographic coordinates for this sphere. 
It is somewhat easier at hrst, however, to parametrize a point on this sphere by its direction 
n, so that functions on phase space are functions of fi. The Q, Weyl, and P symbols for the 
Hamiltonian are 

= x: ^ (A25) 

1=0 m=—i 
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for these in Ref. 
definition 


where are the corresponding symbols for yirn{^)- We gave complete expressions 

29l. but here we only need the asymptotic forms as j ^ oo. Recalling the 

J = J + h (A26) 


we have 


^ /(1 + 0(J 2))W^(h), (A27) 

d)«’^(h) ^ T + 0(r')) Wrn(h). (A28) 

We now observe that + l)Yim{'n) = £^l£m(h), where C = — i(h x Vn) is the angular 
momentum operator (on phase space, and not the quantum mechanical Hilbert space). 
Hence, we may write 

C'’= (lT^ + 0(rt)<. (A29) 

It follows that 

= (l T ^ + 0{r^)')H^, (A30) 

a result which makes no reference to Y^m(h) and is therefore valid independent of the form 
in which the Weyl symbol is given. It has a pleasing similarity to Eq. flAl5l) etc. if we recall 
that —C^ is (the angular part of) the Laplacian on the sphere. 
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